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If the value of On is within the reach of the upcoming generation of long-baseline experiments, 
T2K and NO^A, we show that a low-energy neutrino factory, with peak energy in the few GeV 
range, would provide a sensitive tool to explore CP-violation and the neutrino mass hierarchy. We 
consider baselines with typical length 1000-1500 km. The unique performance of the low energy 
neutrino factory is due to the rich neutrino oscillation pattern at energies between 1 and 4 GeV at 
baselines 0(1000) km. We perform both a semi-analytical study of the sensitivities and a numerical 
analysis to explore how well this setup can measure #13, CP- violation, and determine the type 
of mass hierarchy and the 823 quadrant. A low energy neutrino factory provides a powerful tool 
to resolve ambiguities and make precise parameter determinations, for both large and fairly small 
values of the mixing parameter #13. 



PACS numbers: 14.60.Pq 



I. INTRODUCTION 



o 

o 

During the last several years our understanding of the physics of neutrinos has made remarkable progress. The 
experiments with solar [l|, H, 0, 0, d, @ , atmospheric Q , reactor [|[ , and also long-baseline accelerator P neutrinos 
|' have provided compelling evidence for the existence of neutrino oscillations, implying non zero neutrino masses. The 
present data 1 require two large (612 and #23 ) and one small (#13) angles in the lepton mixing matrix and at least 
two mass squared differences, Am, 2 ; = rnj — mf (where rn^'s are the neutrino masses), one driving the atmospheric 
(Amjj) and the other one the solar (Am 2 ,,) neutrino oscillations. The mixing angles 612 and #23 control the solar 
and the dominant atmospheric neutrino oscillations, while #13 is the angle constrained by the data from the CHOOZ 
. i-H , and Palo Verde reactor experiments [H,J3 ■ 

The Super-Kamiokande [3] and K2K [^1 d ata are well described in terms of dominant — > v T (v^ — » v T ) vacuum 
oscillations. The MINOS Collaboration [l5| has recently reported the first neutrino oscillation results from 1.27 x 
10 20 The value of the oscillation parameters from MINOS are consistent with the ones from K2K, as well as 

from SK data. A recent global fit [13] (see also Ref. [HI) provides the following 3cr allowed ranges for the atmospheric 
mixing parameters 

lAmjjjJ = (1.9-3.2) x 10~ 3 eV 2 , 0.34 < sin 2 9 23 < 0.68 . (1.1) 

The sign of Am 2 1; sign(Am 2 1 ), cannot be determined with the existing data. The two possibilities, Am^ > or 
Am 2 ^ < 0, correspond to two different types of neutrino mass ordering: normal hierarchy and inverted hierarchy. In 
addition, information on the octant in which #23 lies, if sin 2 2#23 7^ 1, is beyond the reach of present experiments. 

The 2-neutrino oscillation analysis of the solar neutrino data, including the results from the complete salt phase of 
the Sudbury Neutrino Observatory (SNO) experiment 0, in combination with the KamLAND spectrum data [lj|, 
shows that the solar neutrino oscillation parameters lie in the low-LMA (Large Mixing Angle) region, with best fit 
values Am^ = 7.9 x 10~ 5 eV 2 and sin 2 6 12 = 0.30. 



1 We restrict ourselves to a three-family neutrino analysis. The unconfirmed LSND signal |ldl| cannot be explained within this context 
and might require additional light sterile neutrinos or more exotic explanations. The ongoing MiniBooNE experiment 11] is going to 
test the oscillation explanation of the LSND result. 



2 



A combined 3-neutrino oscillation analysis of the solar, atmospheric, reactor and long-baseline neutrino data [17j | 
constrains the third mixing angle to be sin 2 (9 13 < 0.041 at the 3er C.L. However, the bound on sin 2 # 13 is dependent 
on the precise value of lAmjjjJ, being stronger for larger values of |Am 31 |. The future goals for the study of neutrino 
properties in neutrino oscillation experiments is to precisely determine the already measured oscillation parameters 
and to obtain information on the unknown ones: namely #13, the CP- violating phase 6 and the type of neutrino mass 
hierarchy (or equivalently sign(Am§ 1 )). 

In the next sections we will explore in detail the possible measurement of the two unknown parameters #13 and 5 
with a future neutrino factory [20l| fa cility as thi s appea rs to be amon g th e most promising ways to unveil neutrino 
mixing and leptonic CP violation [H [|| H, H SI S M, M, El, M, M M, H M, EH 2 . 

A neutrino factory consists of a high intensity muon source 3 , an acceleration system, and a muon storage ring 
with long straight sections. Muons decaying along the straight sections create high intensity neutrino beams which 
have a precisely-known flux, divergence, energy spectrum, and neutrino flavor content. The flavor composition of 
the beam depends on whether positive or negative muons are stored in the ring. Suppose, for example, that positive 
charged muons have been stored. Muons decaying in the straight sections will produce beams containing 50% muon- 
antincutrinos and 50% electron- neutrinos: /i + — > e + + i> e + P M . Charged current interactions of the in a distant 
detector will produce /i + ("right-sign" muons, i.e. muons of the same charge as those stored in the neutrino factory). 
In contrast, if the v e 's oscillate to v^s and then interact in the far detector they will produce \T~ ("wrong-sign 
muons" [20I I2H]). Thus, wrong-sign muons are an unambiguous proof of electron neutrino oscillation, in the v e — ► 
channel. This has been called the "golden channel" [26[, and is central to the present study. A magnetized detector 
with excellent charge identification is necessary to exploit the golden channel. 

In the following we will consider a low energy neutrino factory where the stored muons have an energy of 4.12 GeV. 
This is motivated by recent progress in developing a viable concept for a neutrino factory detector with a threshold 
for reliably measuring the sign of muons with momenta down to a few hundred MeV/c |48| . We explore the impact of 
analyzing the "wrong-sign" and "right-sign" muon rates for several energy bins, and consider two reference baselines: 
1280 Km, the distance from Fermilab to Homestake, and 1480 Km, the distance from Fermilab to Henderson mine. 
Our results can be easily generalized to other baselines in the 1200-1500 km range. We find that a simultaneous fit to 
the energy-dependent rates provides a powerful tool to resolve ambiguities and make precise parameter determinations, 
for both large and fairly small values of the mixing parameter #13. 

II. FORMALISM 

In the present study we focus on the capabilities of a low-energy neutrino factory, if the value of #i 3 is within reach 
of the upcoming generation of long-baseline experiments, T2K and NOi/A. 

For neutrino energies E <; 1 GeV, #13 within the present bounds [Tt], EH, and baselines L ^ 0(1000) km [5l| . 

the oscillation probability P (L) can be ex pan ded in the small parameters #13, A12/A13, A12/A and Ai 2 Z/ , where 
A12 = Aml 1 /(2E) and A i3 = Am| 1 /(2S) [26| (see also Ref. [z^): 

{ P\L) * sin 2 23 sin 2 20 13 (^^sm 2 

+ cos 6 13 sin 20 13 sin 20 23 sin 20 12 ^f^ffe sin (^) sin ( (A ^ 13>I ) cos (^f± T 8) ( 2A ) 
+ cos 2 23 sin 2 20 12 (Ai a ) 2 s i n 2(AL) ) 

where the first, second and third terms have been dubbed atmospheric, interference and solar terms, respectively. 
In the following analytical study, we use the constant density approximation for the index of refraction in matter 
A = V2G F n e (L). Here, n e {L) = 1/L / Q n e (L')dL' is the average electron number density, with n e {L) the electron 
number density along the baseline. 

In order to study the sensitivity to CP-violation, we introduce the weighted probability difference between the case 



2 For the prospects of future measurements of these two oscillation parameters at (3 beam experiments I36t l37ll , see Refs. [38l, l39t Bol. |41|. 

mill- 

3 A neutrino factory muon source is an attractive stepping-stone towards a high energy muon collider. 
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of 6 ^ and the one with no CP-violation (5 = 0) 4 : 

(P(L,S)-P(L,0)) 2 
S(S) = V __ 1 , (2 . 2) 

*W = * pcmT^ ■ (2 - 3) 

The quantity S(S) (S(8)) is useful to get an estimate of the energy range, for a fixed baseline, for which the sensitivity 
is maximal. Using Eq. (|2.ip , we find that: 



4 cos 2 6»i 3 sm z 26»i2 cos^ 6*23 -r— — -- cos(<S — ) - cos — — 

S(S)(S(6)) = VAl3/ V 2 / , ± 2 —— (2.4) 



l + 2ecos(<5 ^-) + e 

where we have approximated cos #13 ~ 1. The quantity e is defined as: 

_ cos 6*23 sin 26»i2 A i2 A13L A =p A i3 1 



sin6»23 sin26li3 A13 2 A i3 sin((A T A ls )L/2) : 



(2.5) 



where =F refers to neutrinos (antineutrinos) respectively. At leading order we can neglect ^4/Ai3 terms and e is the 
same for neutrinos and antineutrinos. 

The sensitivity to CP-violation is when the interference term in the oscillation probability, e.g. the second term 
in the r.h.s. of Eq. (|2.ip . cancels. For the values of 6*13 of interest, this happens at the oscillation minima: 

(A13 TA)L/2 = mr n = 1,2,3.... (2.6) 

Neglecting the small correction due to matter effects, we find that the minima in the sensitivity to CP-violation 
corresponds to an energy of E m = 1.4 (1.2) GeV for L=1480 (1280) km, with n = 1. Here and in the following, 
we used Ato^ = 2.4 x 10~ 3 eV 2 . Matter effects modify this result by 10-20 %. At smaller energies additional 
minima and maxima with a fast oscillatory behaviour are present but we will not consider such energy range due 
to detector resolutions, efficiencies and thresholds. As at high energies S(5) and S(6) drop as E~ 2 , we expect a 
maximum in sensitivity at an energy of few GeV. In the case of neutrinos and normal hierarchy (and, for negligible 
A/A13, antineutrinos and inverted hierarchy), an additional minimum in CP-violation sensitivity is found for cos(5 — 
A13L/2) = cos(Ai 3 L/2). This equation has a solution for: 



A13L 6 



2 2 
A 13 i S 



for 5 > 0, (2.7) 



= 7T - - for S < 0. (2.8) 
& 2 

For (5 = 7r/2 (— 7r/2), the energy of the minimum in sensitivity is at E m — 5.7 (1.9) GeV for L = 1480 km and at 
E m = 5.0 (1.7) GeV at L = 1280 km. Let us notice that, if S < 0, the two minima in Eqs. (pT6j) and |2~8)l are 
very close and the maximum in between cannot be fully exploited due to limited energy resolution. The maximum 
in sensitivity at higher energy is located close to the minimum in Eq. (|2.6p . that is in the energy range E~m ^ (2~ 
3) GeV, depending on the value of 8. For non-negative values of 6, it is possible to make full use of the maximum in 
sensitivity between the two minima in Eqs. (|2.6p and (|2.8p . Again, the maximal sensitivity will be achieved at Em 
of few GeV. In both cases, we can obtain a more precise estimate for Em by neglecting e and solving the equation 
(cos{S- AMi)-cos Ayi^ + AM^ sm (£_ Auii) +S i n Aul^ = 0. Notice that this equation holds only far away from 

the minima. Typically, the maximum is reached in the energy range (1.6-5.2) GeV [(1.4-4.5) GeV] for L = 1480 km 
[L = 1280 km]. In particular, for 5 = tt/2 (-tt/2) we have E M ^ 1.7 (3.2) GeV [1.5 (2.6) GeV] with L = 1480 km 
[L = 1280 km]. Matter effects included in e in the denominator modify these results by a 10-20 %. 



4 A similar analysis with similar results could be carried out for the case 8 = it. We present here the analytical study assuming that S is 
within the interval [— tt/2, 7t/2]. Our numerical simulations will consider the full 8 range [— 7T,7r]. 
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For neutrinos and inverted hierarchy, (and antineutrinos and normal hierarchy, neglecting terms of order A/A13), 
we find a similar behaviour. The minimum in addition to the one in Eq. (|2.6[) is reached when cos(<5 + A 13 L/2) = 
cos(Ai 3 L/2) whose solutions are: 

for 5 > 0, (2.9) 
for 8 < 0. (2.10) 

We can also estimate the position of the maximum by solving the equation ^ cos (S + Al 2 3L ) — cos Al 2 L — Al 2 £ ( sin (5 + 

Al 2 £ ) — sin Al £ L ^j = 0. The maximum is reached at approximately the same energy as for neutrinos but with opposite 

d. As far as e is negligible, the sensitivity does not depend on $13. Our approximated analytical results on the energy 
maxima hold for sin 2 #13 10~ 3 (A13L/2)/ sin(Ai3L/2). In conclusion, our analytical study suggests that maximal 
sensitivity to CP-violation is reached in the few GeV range. Notice that, given one type of hierarchy, neutrinos and 
antineutrinos have similar behaviour but for opposite values of S. The combination of the two channels allows to 
reach optimal sensitivity independently of the true value of the CP-violating phase. 

Similar considerations hold also for the sensitivity to the type of hierarchy. We can study a similar quantity: 
(P + — P-) 2 /P + for neutrinos and antineutrinos. We find that the dominant term is proportional to AL while CP- 
violating terms constitute a correction at most of 20%-30% for the highest allowed values of sin 2 2#i3. A minimum in 
the sensitivity is found in correspondence to the minima of the oscillation probability as in the case of CP-violation 
studied above. The sensitivity to the type of hierarchy depends on the value of the 6 phase once the CP-violating 
corrections are taken into account. For < <5 < n/2, the energy for which a maximum in sensitivity is obtained, Em, 
will be an decreasing function of the #13 mixing angle. For example, for 6 = 0, the maximum in sensitivity will be 
reached in the 3.7-2.3 GeV range for sin 2 613 = 0.01-0.1. Conversely, for n/2 < 6 < n, Em will increase with #13 but 
typically remain in the 1.3-1.6 GeV range. We can conclude that maximal sensitivity is reached for energies around 
1.3-4 GeV. For antineutrinos, a similar behaviour can be found, with the exchange of 5 in n — 8. 

In order to study the sensitivity to CP-violation and type of hierarchy, by exploiting the number of events in a 
simulated experiment, the energy dependence of the cross sections and fluxes should be included. We can expect the 
value of the energy for which the maxima of sensitivity are reached to be shifted slightly at higher values. 

From the above considerations, we can conclude that the use of a detector with a low threshold and good energy 
resolution and efficiency at E <; 1 GeV is crucial for exploiting the potentiality of a neutrino factory with baselines in 
the 1000-1500 Km. In addition, a high energy neutrino beam is not necessary and it is sufficient to use lower energies 
with respect to the commonly studied options for a neutrino factory with muon energies of 20 or 50 GeV. 



A 13 L 6 

— =1T -2 
A13L = 5_ 

2 2 



III. THE LOW ENERGY NEUTRINO FACTORY CONCEPT 



A Neutrino Factory consists of (i) an intense low energy muon source, (ii) a muon beam manipulation and cooling 
system to maximize the number of muons within a useful acceptance, (iii) a "pre- accelerator" to accelerate the muons 
from low kinetic energies (typically 100-200 MeV) to about 1 GeV, (iv) a system of one or more accelerators to 
further accelerate the muons to the desired final energy, and (v) a muon storage ring with long straight sections. 
Design studies [3, EH EH have shown that, for a 20 GeV Neutrino Factory, the 1-20 GeV acceleration systems 
are expected to account for about 26% of the estimated cost. Hence, if the physics goals can be met using muons 
with energies much lower than 20 GeV, there is a significant cost advantage. In the following, we first discuss the 
performance of the far detector (which places a lower limit on the desired Neutrino Factory muon energy), and then 
discuss the low energy Neutrino Factory and its performance. The primary neutrino oscillation channel at a Neutrino 
Factory requires identification of wrong-sign muons, and hence a detector with excellent muon charge identification. 
To obtain the required event rates, the far detector fiducial mass needs to be at least O(10 Kt), and therefore we 
require a very large magnetized detector. Early studies [HJ based on a MINOS-like segmented magnetized detector 
suggested that, to reduce the charge mis-identification rate to the 10 -4 level while retaining a reasonable muon 
reconstruction efficiency, the detected muon needs to have a minimum momentum of 5 GeV. The analysis obtained 
a 50% reconstruction efficiency for CC neutrino interactions exceeding ~ 20 GeV. This effectively places a lower 
limit of about 20 GeV on the desired energy of the muons stored in the Neutrino Factory. However, a more recent 
analysis [48| has shown that, with more sophisticated selection criteria, high efficiencies (> 80%) can be obtained for 
neutrino interactions exceeding ~ 10 GeV, with efficiencies dropping to ~ 50% by 5 GeV. The new analysis suggests 
a MINOS-like detector could be used at a Neutrino Factory with energy less than 20 GeV, but probably not less than 
10 GeV. If we wish to consider lower energy Neutrino Factories, we will need a finer grained detector that enables 
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Turn Number 


1 
1 


z 


Q 
O 


A 





Initial Energy (GeV) 


1.0 


1.96 


2.92 


3.88 


4.84 


Final Energy (GeV) 


1.48 


2.44 


3.40 


4.36 


5.32 


f decay = 100m/ 7CT (%) 


1.30 


0.73 


0.51 


0.39 


0.32 


Ndecay per year (xlO 18 ) 


9.8 


5.5 


3.8 


2.9 


2.4 



TABLE I: Useful positive muon decays in one straight section of an RLA designed to accelerate from 1 GeV to 5.8 GeV in 5 
turns. The straight section and arc lengths are, respectively, 100m and 30m. The numbers tabulated correspond to 7.5 x 10 20 
injected muons, or roughly one years operation. 

reliable sign-determination of lower energy muons with good efficiency. One way to achieve this would be to use a 
totally active magnetized segmented detector, for example a NOvA-like detector [47| within a large magnetic volume. 
Initial studies [48| show that, for this technology, the muon reconstruction efficiency is expected to approach unity 
for momenta exceeding ~ 200 MeV/c, with a charge mis-identification lower than 10~ 4 (10~ 3 ) for momenta exceeding 
approximately 400 MeV/c (300 MeV/c). 

Further studies are needed to fully understand the efficiency and mis-identification as a function of neutrino energy, 
but there is hope that a Neutrino Factory detector might be designed with good wrong-sign muon identification and 
high efficiency for neutrino energies as low as 1 GeV, or perhaps a little lower. Given these recent developments, in 
our analysis we will assume that a massive magnetized detector can be designed to identify wrong-sign muons with 
full efficiency for neutrino interactions above 0.8 GeV, and zero efficiency below this energy. We will see that the 
excellent physics capability of a low energy Neutrino Factory motivates striving for a detector that can achieve this 
demanding performance. 

With a magnetized far detector concept that makes plausible precision measurements of neutrino interactions down 
to about 0.8 GeV, we are motivated to consider Neutrino Factories with stored muon energies of a few GeV. In 
present designs for a 20 GeV Neutrino Factory [4q . there are at least two acceleration stages that accelerate the 
muons from about 1 GeV to 20 GeV. Depending on the design, these accelerators consist either of Recirculating 
Linear Accelerators (RLAs) or Fixed Field Alternating Gradient accelerators (FFAGs). A few GeV Neutrino Factory 
would require only one of these acceleration stages. Note that the RLAs have long straight acceleration sections 
which, if pointing in a suitable direction, could provide a neutrino beam with a time-dependent energy that varies 
from 200 MeV up to the final energy. This might facilitate a tunable-energy Neutrino Factory. To illustrate this, 
Table U shows, for 7.5 x 10 20 positive muons (per year) injected into an 1-5.8 GeV RLA, the number of muon decays 
in a given straight section at each intermediate energy. Only 7% of the injected muons decay during the acceleration, 
and hence 93% (7 x 10 20 ) are available to be injected into a dedicated fixed-energy Neutrino Factory. If the Neutrino 
Factory straight section length is 30% of the ring circumference, this would provide an additional 2 x 10 20 useful muon 
decays per year at 5.8 GeV. Note that using the RLA to provide a neutrino beam would provide flexibility in choosing 
the desired neutrino energy spectrum. The acceleration cycle could, in principle, be varied to keep the muons for as 
long as desired at any intermediate energy. Hence, the 2 x 10 20 useful muon decays could be redistributed amongst 
the intermediate energies, as needed. However, the flexibility of using the RLA to provide a low energy neutrino 
beam comes at the cost of a more complicated design for the accelerator. In particular, the angular divergence of 
the beam within the straight section needs to be small compared to the angular spread of the neutrinos generated 
in muon decay. If designs that achieve this prove to be impractical or expensive, flexibility could also be achieved 
by designing the RLA so that the muons could be extracted on any given turn and injected into one of several fixed 
energy Neutrino Factory storage rings (note that the cost of the storage rings is believed to be small compared to the 
cost of the RLA). 

In the following we will show that a low energy Neutrino Factory with a fixed energy of 4.12 GeV would provide a 
sensitive tool for exploring neutrino oscillations if #13 is " large" . This energy would require about 4 turns in a single 
RLA. Note that the sensitivity of a Neutrino Factory experiment depends upon the event statistics, and hence upon 
the product of the detector fiducial mass, the length of the data taking period, and the number of muons per unit 
time decaying in the appropriate Neutrino Factory straight section. Initial studies have considered, as reasonable, 
a totally active magnetized detector with a fiducial mass of about 20 Kt. Present Neutrino Factory studies suggest 
that it would be reasonable to expect, for a Neutrino Factory with (without) a muon cooling channel before the 
pre-accelerator, about 5 x 10 20 (3 x 10 20 ) useful positive muon decays per year and 5 x 10 20 (3 x 10 20 ) useful negative 
muon decays per year in a given Neutrino Factory straight section. Hence, a conservative estimate of the sensitivity 
of a Neutrino Factory experiment might be based upon 5 years data taking with 3 x 10 20 useful muon decays of each 
sign per year in the storage ring, and a detector fiducial mass of 20 Kt, corresponding to 3 x 10 22 Kt-decays for each 
muon sign. A more aggressive estimate might be based upon 10 years data taking with 5 x 10 20 useful muon decays 
of each sign per year in the storage ring, and a detector fiducial mass of 20 Kt, corresponding to 1 x 10 23 Kt-decays 
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for each muon sign. 

IV. NUMERICAL ANALYSIS: DEGENERATE SOLUTIONS 

We can ask ourselves whether it is possible to unambiguously determine S and #13 by measuring the transition 
probabilities v e — > and — > at fixed neutrino energy and at just one neutrino factory baseline. The answer 
is no. It has been shown [2!| that, by exploring the full (allowed) range of the S and #13 parameters, that is, 
— 180° < S < 180° and 0° < 6*13 < 10°, one finds, at fixed neutrino energy and at fixed baseline, the existence of 
degenerate solutions (# 13 , S ), that have been labelled in the literature intrinsic degeneracies, which give the same 
oscillation probabilities as the set (#13, S) chosen by nature. More explicitly, if (#13, S) are the correct values, the 
conditions 

JW*i3> 5 ')= i W*l3,*) 

can be generically satisfied by another set (9 13 , S ). 

Additional solutions might appear from unresolved degeneracies in two other oscillation parameters: 

1. At the time of the future neutrino factory, the sign of the atmospheric mass difference Am^ may remain 
unknown, that is, we would not know if the hierarchy of the neutrino mass spectrum is normal or inverted. In 
this particular case, P(9 13 ,5 , — AmjJ = P(9is, S, Am 31 ). 

2. Disappearance experiments only give us information on sin 2 2023^ is #23 in the first octant, or is it in the second 
one, (n/2 — #23)? • In terms of the probabilities, P(9i 3 , 6 , § — #23) = P(9\3,S, #23)- 

This problem is known as the problem of degeneracies in the neutrino parameter space [l^, [3(J IH, H3, HH HH . All 
these ambiguities complicate the experimental determination of S and #13. Many strategies have been advocated to 
resolve this issue which in general involve another detector [29|, [3l|, HH, [H, [H, HE E3, El| or the combination with 
another experiment [H, HEM HE EE M, EE EI EE EE E3 . 

In the present study we show that, if the value of #13 turns out to be not very small, a low energy neutrino 
factory provides the ideal scenario for the extraction of the unknown parameters as well as for resolving the discrete 
degeneracies. The reason is simple: at distances of 0(1000) km the neutrino oscillation pattern is extremely rich at 
neutrino energies below 4 GeV. We have thus explored a single decaying muon energy scenario = 4.12 GeV. 

By exploiting the energy dependence of the signal, it is possible to disentangle $13 and 5 and to eliminate the 
additional solutions arising from the discrete ambiguities. We have divided the signal in four energy bins. The energy 
binning of the signal has been chosen accordingly to the analytical study. The energy range of these four energy bins 
is [0, 0.8], [0.8, 1.5], [1.5, 3.5] and [3.5, 4.12] GeV. We will show the physics potential of the chosen energy binning here 
in the next subsection. The detection efficiencies are considered as perfect (100%) above the first bin, as described in 
the previous section. 

Two possible baselines have been carefully explored: 1280 Km, the distance from Fermilab to Homestake, and 
1480 Km, the distance from Fermilab to Henderson mine. The results are presented for the two possible scenarios 
described in the previous section, in order to quantify the benefits of increased exposure times and muon intensities: 
a conservative scenario of 3 x 10 22 Kton-decays and a more ambitious one with 1 x 10 23 Kton-decays. 

All numerical results simulated in the next subsections here have been obtained with the exact formulae for the 
oscillation probabilities. Unless specified otherwise, we take the following central values for the remaining oscillation 
parameters: sin 2 9 12 = 0.29, Am^ = 8 x 10~ 5 eV 2 , [Am^l = 2.5 x 10~ 3 eV 2 and 6*23 = 40° . The X 2 for a fixed 
baseline A is defined as: 

X* = E E Hp - N tv) C iM Hp' ~ N L') > 

i,j P,p' 

where is the predicted number of muons for a certain oscillation hypothesis, n^ p are the simulated "data" from 
a Gaussian or Poisson smearing and C is the 2Nbi n x 2Ni,i n covariance matrix, that will contain statistical errors and 
a 2% overall systematic error. All the contour plots presented in the following in a two parameter space have been 
performed assuming 2 d.o.f statistics. 
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A. Optimizing the energy binning 

In this subsection we provide an explanation for the energy binning chosen here, crucial to resolve the additional 
solutions (degeneracies). As first noticed in Ref. [59], the location of the degeneracies is E,L dependent. For large 
#13, the location of the intrinsic degeneracies is given by: 

8 ~ 7T — 8, 

4 ~ 13 + cos5sin20 12 ^^cot0 23 cot(^^) . (4.2) 

Notice that the shift 9 13 — 813 depends on the energy and the baseline through the function cot ^ A "j, 1 ^ ■ If the 

function cot ^— -j-ff— ^ changes sign from one energy bin to another, the degenerate solutions will appear in different 

regions of the parameter space. Consequently, the combination of fits to the various energy bins can eliminate the 
intrinsic degeneracies. 

To illustrate this, Figure ©(a) depicts results from fitting the simulated data for each of the energy bins. The 
simulation is for L = 1480 km, 6*13 = 3°, S = 0°, normal mass hierarchy, and 823 = 40°. The fit results shown in the 
figure correspond to the correct hierarchy and 823 (the impact of including the discrete degeneracies will be discussed 
later). The 90%, 95% and 99% CL contours resulting from the fits are shown for the second energy bin (blue), third 
energy bin (cyan,) and fourth energy bin(magenta). Notice that in addition to the correct solution, there are also 
fake solutions for which 8 ~ ir, as indicated by Eq. (]4 . 2[) . However, the fake solutions from the fit to the third energy 
bin get opposite displacements 13 — 613 from those from the fits to the second and fourth energy bins. The relative 
displacement (positive or negative) is given by the sign of the trigonometric function cot ^ A 4 ™J^ L J , where (E) is the 
medium energy for a fixed bin. A combination of fits to the second, third and fourth energy bins will therefore help 
in resolving the intrinsic degeneracies. 




FIG. 1: 90%, 95% and 99% (2 d.o.f) contours resulting from fits to simulated data for L = 1480 km, #13 = 3°, and 8 = 0° 
(position in parameter space denoted by a star). The blue, cyan and magenta contours represent the results from fits to the 
second, third and fourth energy bins, respectively. In the left panel, we have neglected the impact of the discrete degeneracies. 
In the right panel, the fit assumes incorrectly that the sign of the atmospheric mass difference is negative. The atmospheric 
mixing angle is fixed to #23 = 40° . 



A similar exercise to the one described above can be done for the discrete degeneracies. For instance, the wrong- 
sign(Am| 1 ) additional solutions will be located in different regions of the parameter space for each energy bin, see 
Fig. (HJ)(b), and therefore fits to the combination of the energy bins will result in a resolution of the sign degeneracies. 
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The location of the fake solutions for the second energy bin is different from those for the third and fourth energy 
bins, since at lower energies matter effects are less important. 

Resolving the additional #23-octant degeneracy is, in general, very difficult. As shown in Ref. [5§|, for large #13, the 
location of the #23 degeneracies is given by: 

sin S ~ cot 823 sin 5, 

sin2^ 2 ^M 1 / / AmU\ ( , Am 2 L\\ , . 

tfx, - tan e 23 d 13 + ^a^j ( cos ( s - -dr) - tan * 23 cos (* --uh))- (43) 

This system describes two solutions. For one of them, the L and E dependent terms in Eq. (|4.3[) tend to cancel 
for 023 ~> t/4, resulting in # 13 = #13 and <5 = <5 in this limit. The second solution coincides in this limit with the 
intrinsic degeneracy, Eq. (|4.2p . Notice that no fake solutions are expected for | cot #23 sin <5| > 1. Figure @ illustrates 
the equivalent exercise to those performed above for the intrinsic and for the wrong-sig^Am 2 -^ degeneracies. The 
simulated data is for the #23 in the first octant, i.e. 6*23 = 40°, while the fit is performed incorrectly assuming the second 
octant, i.e. #23 = 50°. Notice from the results depicted in Fig. ©(a) that there are two sets of degenerate solutions, 
as previously discussed: those which resemble the correct values and those which are related to the intrinsic solution. 
While the location of the former is E, L independent, the location of the latter will depend on E, L, and therefore the 
combined fits to the various energy bins will help in resolving these degeneracies. Note that the degeneracies which 
are closer to the correct values are extremely difficult to resolve. The information in the second bin is crucial: if the 
detector efficiency in the second bin ([0.8, 1.5] GeV) is sufficiently high (we are assuming 100%), the combination of 
fits to the various energy bins will resolve the additional solutions related to the wrong choice of the atmospheric 
mixing angle octant. Since the second bin is the lower energy bin, it turns out that for these lower energies the solar 
term, see Eq. (I2.1|) . is the dominant one for #13 < 3°, larger than both the atmospheric and the interference terms. 
The solar term goes as cos 2 6*23, while the atmospheric term goes as sin 2 #23, therefore, exploiting this low energy 
bin is crucial to resolve the 6*23 degeneracy. The resulting fit, after the combination of the data in the three energy 
bins considered here, is degeneracy free down to very small values of #13 ~ 1° (see Fig. ©(b)), since the additional 
solutions from fits to the second bin lie on different locations in parameter space than those for the third and fourth 
energy bins. 




FIG. 2: (a) 90%, 95%and 99% (2 d.o.f) contours resulting from the fits to the data at L — 1480 km, assuming that the nature 
solution is 813 = 3° and 5 = 0° (denoted by a start). The blue, cyan and magenta contours represent the results with the second, 
third and fourth energy bin data, respectively. We assume that the true value of the atmospheric mixing angle is 823 = 40° 
(first octant) but the data is fitted to 823 = 50° (second octant), (b) Same than in (a) but assuming that the nature solution is 
813 = 1° and 8 — 0° (denoted by a start). 
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B. Exploring the disappearance channel 

We explore here the measurement of the atmospheric neutrino oscillation parameters, 6*23 and Ato^ making use of 
the disappearance channel. We study as well the impact of this channel regarding the # 2 3-octant degeneracy. The 
disappearance channel at the neutrino factory has been already considered in the literature [24], HH and it has been 
widely and carefully explored in Ref. {69} . The vacuum — > oscillation probability expanded to the second order 
in the small parameters #13 and (A 12 L/£ ; ) reads [7(| 

Ptyn ~+v») = 1 - [sin 2 26*23 - s 2 3 sin 2 26 13 cos26> 23 ] sin 2 1 A ' J 



2 



A 12 L 



[s\ 2 sin 2 26*23 + Jsl 3 cos 5] sin(A 23 .L) 



2 

A L\ 2 

J ' * [4, sin 2 20i2 + s 2 2 sin 2 26» 23 cos(A 23 i)] , (4.4) 

where J = 0036*13 sin 2#i2 sin 26*13 sin 26*23 and A23 = Am 32 /2E, A12 = A?7i|i/2.E. The first term in the first 
parenthesis is the dominant one and is symmetric under 6*23 — > n/2 — #23- However, when a rather large non- vanishing 
6*i3 is switched on, a 6'23-asymmetry appears in Eq. (|4.4[) and the octant in which 6*23 lies can be extracted from 
disappearance data, as will be shown in our numerical results. We assume here the same detection efficiencies 5 and 
energy binning than those which will be considered for the golden v e — > transition. A global 2% systematic error 
has been included in the \ 2 fits to the atmospheric neutrino parameters. We present our results in the (sin 2 #23, Am 2 ^) 
plane in Figs. (O for two simulated values of 6*13, and two simulated values for sin 2 6*23: sin 2 #23 = 0.4 and sin 2 6*23 = 
0.44. The detector is located at the Henderson mine at a baseline of L = 1480 km (similar results are obtained for 
the Homestake detector location). The CP violating phase S has been set to zero. Notice that this channel is able 
to reduce atmospheric parameter uncertainties to an unprecedented level: the resolution in sin 2 6*23 is astonishing, 
maximal mixing can be excluded at 99% CL if sin 2 6*23 < 0.48 (0 23 < 43.8°), independently of the value of 6>i 3 . In 
addition, for a relatively large value of 6*13 > 8°, the 6 ) 2 3-octant degeneracy will not be present at the 99% CL for 
sin 2 6*23 < 0.44 (6^3 < 41.5°), if #13 is treated as a fixed parameter. These results have been obtained for the more 
conservative neutrino factory scenario described above, a scenario with 3 x 10 22 Kton-decays for each muon sign. 
Since the statistics and the size of the expected signal are both large in disappearance measurements, the error on 
the parameters will be dominated by the systematic error and a more ambitious scenario with higher statistics (with 
1 x 10 24 Kton-decays) will not improve much these results. 



C. Simultaneous fits to #13 and S 



In this subsection we exploit the golden channel, i.e. the v e (y e ) — > ^(z/^) transitions to extract the unknown 
parameters 6*13 and S. We start exploring the more conservative neutrino factory scenario with 3 x 10 22 Kton-decays. 
We present in Figs. ((4]) the 90%, 95% and 99% CL contours for a fit to the simulated data from a future low energy 
neutrino factory with the detector located at Homestake, at a baseline L = 1280 km (left panel) and at Henderson, 
at a baseline L — 1480 km (right panel). The "true" parameter values that we have chosen for these examples are 
depicted in the figures with a star: we have explored four different values of S = 0°, 90°, —90° and 180° and #13 = 8°. 
The simulations are for the normal mass hierarchy and #23 in the first octant (sin 2 6*23 = 0.41 which corresponds to 
6*23 = 40°). Our analysis includes the study of the discrete degeneracies. That is, we have fitted the data assuming 
both the wrong hierarchy and the wrong choice for the #23 octant (i.e. negative hierarchy and sin 2 #23 = 0.59, which 
corresponds to 6*23 = 50°) and the additional solutions (if present) will be shown in red and in cyan, respectively. 
Notice that in Figs. ([4} the sign ambiguity is solved at the 99% CL. The additional solutions associated to the wrong 
choice of the 6*23 octant are not present at the same CL due to the information extracted from the disappearance 
channel. 

Similar results are obtained for smaller values of 6*13, see Figs. ([5]). Notice that the performance of the low energy 
neutrino factory is unique: the sign(Am 2 i) can be determined at the 99% if 6*13 > 2° independent of the value of the 
CP phase S. Regarding the 6 ) 2 3-octant ambiguity, it can be removed at the 99% CL down roughly to #13 > 1° for a 



We believe this is conservative since less aggressive cuts are required to reduce backgrounds for the disappearance channel than those 
required for the appearance channel. 
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sin 8 23 sin a 23 

(a) 0i 3 = 8° (b) 6>i3 = 4° 

FIG. 3: 90%, 95% and 99% (2 d.o.f) CL contours resulting from the fits at L — 1480 km assuming two central values for 
sin 2 023 = 0.4 and 0.44 and Am^ = 2.5 x 10~ 3 eV 2 . In the right (left) panel, 813 = 8° (4°). A larger value of 813 introduces a 
larger asymmetry and the four-fold degeneracy in the atmospheric neutrino parameters is solved. The statistics considered for 
both simulations corresponds to 3 x 10 22 Kton-decays. Only disappearance data have been used to perform these plots. 

nature's choice of sin 2 6*23 = 0.41, independent of the value of S, provided that the conservative estimate of 3 x 10 22 Kt- 
decays for each muon sign are feasible. The 623 octant degeneracy is solved with the information contained in the 
second energy bin data, which is sensitive to the solar term, as mentioned before. 

Since the results are very similar for the two baselines explored here, the physics reach with the more aggressive 
estimate of 1 x 10 23 Kt-decays for each muon sign is illustrated for only one baseline, L — 1480 km (Henderson mine 
site) and for smaller values of 6*13. Figure ([6]) shows fit results for two simulated values of #13 (2° and 1°). The mass 
hierarchy can be determined at the 99% CL if 813 > 1° independent of the value of the CP phase S. In addition, for 
our example with sin 2 6*23 = 0.41, the 6*23 octant ambiguity can be resolved at 99%CL for all values of the CP phase 
S provided 6 13 > 0.6°. 

We summarize the reach of the low energy neutrino factory with two exclusion plots which illustrate the performance 
of the experiment explored here. We have taken into account the impact of both the intrinsic and discrete degeneracies 
to depict the excluded regions. For both exclusion plots we have assumed that the detector is located at the Henderson 
mine at L = 1480 km. The results for the closer baseline (L — 1280 km) are very similar. 

Figure (O depicts the region in the sin 2 26*13 , fraction of S plane for which the hierarchy can be resolved at the 95% 
CL assuming 2 d.o.f statistics, for both scenarios, the more conservative one, in which the exposure is 3 x 10 22 Kton- 
decays, and the more aggressive scenario in which the exposure is 1 x 10 23 Kton-decays. Notice that the hierarchy 
could be determined in both scenarios if sin 2 26*13 > 0.01 (i.e. 6*13 > 3°) regardless of the value of the CP violating 
phase 5. 

Figure (O depicts the region in the (sin 2 26>i3, 8) plane for which a given (non-zero) value of the CP violating phase 
can be distinguished at the 95% CL from the CP conserving case, i.e. S = 0, ±180° (assuming 2 d.o.f statistics). This 
exercise is illustrated for the two scenarios considered in this study, the more conservative one, in which the exposure 
is 3 x 10 22 Kton-decays, and the more agressive scenario in which the exposure is 1 x 10 23 Kton-decays. Notice that 
the CP violating phase 5 could be measured with a 95% CL error lower than 20° in both scenarios if sin 2 26*13 > 0.01 
(i.e. #13 > 3°), reaching an unprecedent precision for larger values of #13. For smaller values, 0.001 < sin 2 26>i3 < 0.01, 
and in the more conservative scenario, the presence of the sign-degeneracy compromises the extraction of the CP 
violating phase 5. On the other hand, in the more aggressive scenario, a CP violating effect could be established at 
the 95% CL if sin 2 26 13 > 0.001 for 20° < S < 160° (-160° < 8 < -20°). 
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(a) Homestake, L = 1280 km. 




(b) Henderson, L = 1480 km 



FIG. 4: 90%, 95% and 99% (2 d.o.f) CL contours resulting from the fits at L = 1280 km (left panel) and L = 1480 km (right 
panel) assuming four central values for S = 0° , 90° , —90° and 180° and #13 = 8°. The statistics considered for both simulations 
corresponds to 3 x 10 22 Kton- decays. 




(a) Homestake, L = 1280 km. 




(b) Henderson, L = 1480 km 



FIG. 5: The same as Figs. but for 613 — 3°. The additional solution (only at the 99% CL) associated to the wrong choice 
of the neutrino mass ordering is depicted in red. The additional solutions (only at the 99% CL) arising from the wrong choice 
of the 823 octant are depicted in cyan. The statistics considered for both simulations corresponds to 3 x 10 22 Kton-decays. 



V. CONCLUSIONS 



We have shown here the enormous physics reach of a novel neutrino factory concept, a low energy neutrino factory, 
in which the stored muons have an energy of 4.12 GeV. 

We have exploited both the disappearance (v^ — ► v^) and the golden (y e — > v^) channels by measuring the "right- 
sign" and the "wrong-sign" muons at a two possible baselines: 1280 Km, the distance from Fermilab to Homestake, 
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FIG. 6: 90%, 95% and 99% (2 d.o.f) CL contours resulting from the fits at L = 1480 km assuming four central values for 
8 — 0°, 90°, —90° and 180° and 8n = 2° in the left panel (813 = 1° in the right panel). The additional solutions associated to 
the wrong choice of the neutrino mass ordering are depicted in red. The statistics considered for both simulations corresponds 
to 1 x 10 23 Kton-decays. 



and 1480 Km, the distance from Fermilab to Henderson mine. The results presented here can be easily generalized 
to other baselines in the 1200-1500 km range. 

We illustrate the results of the analysis of the energy binned signal for a facility with (a) 3 x 10 22 Kt-decays for 
each muon sign and (b) 1 x 10 23 Kt-decays for each muon sign. The novel setup presented here could extract the #13 
angle, the neutrino mass hierarchy and the leptonic CP violating phase 8 with unprecedented precision. 

The unique performance of the low energy neutrino factory (when compared to the common 20 — 50 GeV neutrino 
factory) is due to the rich neutrino oscillation pattern at energies between 1 and 4 GeV at baselines 0(1000) km. Recent 
studies have shown that it could be possible a Neutrino Factory detector with good wrong-sign muon identification 
and high efficiency for neutrino energies as low as 1 GeV, or perhaps a little lower. Therefore, to evaluate the physics 
potential of a low energy neutrino factory, we have assumed 100% efficiency above a threshold energy of 0.8 GeV, 
and zero efficiency below this threshold. This naive model for the detector performance will need to be updated once 
further work has been done to better understand the expected detector energy dependent efficiency. 

With this caveat we find that maximal atmospheric neutrino mixing can be excluded at 99% CL if sin 2 23 < 0.48 
(623 < 43.8°). If the atmospheric mixing angle is not maximal, for a nature's choice of sin 2 6*23 = 0.4, the octant in 
which 023 lies could be extracted at the 99% CL if 13 > 1° (0 13 > 0.6°) with an exposure of 3 x 10 22 Kt-decays 
(1 x 10 23 Kt-decays) for each muon sign, independently of the value of the CP violating phase 8. The neutrino mass 
hierarchy could be determined at the 95% CL, and the CP violating phase 8 could be measured with a 95% CL 
error lower than 20°, if sin 2 26*13 > 0.01 (i.e. 6*13 > 3°) assuming the more conservative exposure scenario. All the 
sensitivities quoted here are computed assuming the 2 d.o.f statistical approach, and in our analysis we have included 
statistical and a 2% overall systematic error. 

In summary, the low energy neutrino factory scenario could provide the ideal laboratory for precision lepton physics 
if the mixing angle #i 3 > 2°. 
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